Spiky Phases of Smooth Membranes. 
Imphcations for Smooth Strings 
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We point out a possible mechanism by which smooth surfaces can become spiky as the constant 
of curvature stiffness k falls below certain critical values. This happens either in a single first-order 
transition, or in a sequence of two Kosterlitz-Thouless-like transitions. There may also be additional 
phases in which the spikes form a hexagonal solid-like array or a disordered liquid-like structure. 
Our discussion suggests that there exist smooth strings between quarks. 
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I. INTRODUCTION 



We would like to point out a possible rather universal mechanism by which spiky superstructures can form on a 
membrane. These can undergo several interesting phase transitions, whose nature we predict. 

The existence of superstructures on membranes has been investigated by many authors, in particular wormhole |^ 
and egg-carton shapes |^,||. There are also experimental indications that such structures exist (see Fig. [^). 




FIG. 1. Rough Surface Structure observed by B. Kloesgen and W. Helfrich [g] in dioleoylphosphatidylcholine (DOPC) 
bilayers. 

Consider closed amphiphilic vesicles dispersed in water, forming in general smooth tensionless surfaces |^] whose 
effective energy is governed by curvature stiffness S: 



^ J d^^Vg ho(ci + C2)^ + K0C1C2] , 



(1) 



where kq is of the order of eV. Here (ci + C2)/2 and C1C2 denote the mean and the Gaussian curvature, respectively, 
X = (cc^jX^) are arbitrary parameters of the surface, and gij{x) is the intrinsic metric. This energy is certainly only 
an approximation, valid for small curvatures. If the curvatures increase, there will be deviations from this simple 
quadratic behavior. Moreover, since membranes are composed of rod- like molecules of a certain length tq, there 
exists a natural maximal curvature Ci = C2 = 1 /rQ beyond which a membrane cannot be bent without destroying the 
microscopic structure. Between zero and this maximal value, the increase in energy will slow down. This property 
is caused by various important contributions to the bending stiffness |^,^. As a function of the mean curvature, the 
bending energy of many membranes may have a maximum at a certain value Cm of (ci -|- C2)/2, as sketched in Fig. |^. 
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(ci +C2)ro/2 

FIG. 2. Possible behaviors of bending energy as a function of mean curvature (ci + C2)/2 measured in units of the inverse 
molecular size ro- 

Such a maximum supplies the system with a second characteristic length Tm = l/cm at which interesting new 
phenomena should be observed. A further length scale arises from higher powers of the Gaussian curvature in the 
energy density. Assuming, for example, the presence of terms K2(ciC2)^ + ^4(0102)'*, with coefficients K2 ~ — eVA^, 
K4 «eV (30A)®, Monte-Carlo simulations |9| have shown the existence of a periodic egg carton superstructure with a 
period of the order of = 6OA. 

Membranes with an additional length scale of either type will be the objects of study in this note. Certainly, the 
energy of a single spike will be quite large (maybe of the order of lOfc^T) so that their Boltzmann factor is quite small. 
We shall see, however, that configuration entropy can compensate this large energy at high enough temperatures. 



II. REMINDER OF FLUCTUATION PROPERTIES OF MEMBRANES 



2. Before we come to our actual theory, let us first recall some well-known properties of membrane fiuctuations 
governed only by the energy (|l|), with a partition function (in natural units where the Boltzmann constant is equal 
to 1) 



Z^^e-^/^, (2) 

conf 



where Y^o^f denotes a sum over all membrane configurations. As a consequence of the anharmonic nature of the 
energy (n^, when expressed in terms of the positions of the membrane molecules, thermal fiuctuations make the first 
coupling constant, the extrinsic curvature stiffness, soften with temperature as ||lO| , |rT|| 

« = Ko-3(r/87r)log(riVr2^), (3) 

Here ryv is a short-distance (UV) cutoff equal to the size of the molecules tq, whereas hr, is a long-distance (IR) cutoff 
set by the size of a vesicle. The small parameter in the perturbation expansion leading to (||) is the inverse stiffness 
the flexibility. The result (^ is obtained from an infinite bubble sum of diagrams, which is equivalent to a 
Hartree-Fock-Bogoljubov approximation [pl]| . That approximation gives usually a good idea also for large fiexibilities, 
i.e. small k. If we trust Eq. (^ in this regime, we conclude that the extrinsic curvature vanishes for large vesicles, 
whose size exceed the so-called persistence length [ p^ 

e-ruve^-'-''/^^. (4) 

On the basis of this, we may expect vesicles of a size much larger than ^ to look crumpled. This has indeed been 
confirmed in computer simulations 

For the second coupling constant kq in (|^), the Gaussian curvature stiffness, this effect is absent since this constant 
hardens as follows |13| 

At = So + (10/3)(T/8^) log {rjjrl^) . (5) 

An important question is whether the persistence length (|^) refiects a true property of the theory, or is merely a 
consequence of the one- loop approximation which will not survive higher- loop corrections. In an attempt to answer 
this question, consider a single surface fluctuating around an infinitely large planar configuration in d dimensions with 
periodic boundary conditions. This surface may be described by a vector field X^^{x) {d = 1, . . . ,d), giving rise to 
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the intrinsic metric gij{x) = diX^{x)djX^{x) where = 1,2). The contraction of the spatial vector indices ^ is 
performed via the euclidean unit matrix G^j^ = Sf^^- In terms of X^{x), the energy (|l]) can be written more exphcitly 
as 

E=\j <fx^^) [KoD'^Xt'{x)D'^X^,{x) + Ko {d''X^D''X^,-D,DxX^'D''D^X^,)] , (6) 

where I?^ is the covariant derivative formed with the help of the Christoffel symbol Tij'^ = g^^{digji +djgii — digij)/2. 
The Gaussian curvature energy can be ignored, since it is a constant depending only on the genus of the surface. 

We shall choose a special parametrization due to Gauss, in which X^{x) = x^, X'^{x) = x^, and X°+^(a:) = u°'{x) 
with (a = 1, . . . , d — 2), so that gij{x) = Sij + diUa{x)djU°' (x) , with the contraction over indices a being performed 
via the (d — 2) x (d — 2)-dimensional submatrix Sat of G^^. Then we rewrite the energy (^) in yet another form 

e' = yJ '''^^ [D^x'^D^Xa + y'{d^x"-d,Xa + % - K,)] , (7) 

where hij{x) and Xij{x) are two auxiliary symmetric fields. The partition function of the surface is then given by the 
functional integral 

Z = I VuVX^jVhije^^'/'^. (8) 



The integral over Xij (x) ensures that the auxiliary field hij (x) coincides with the induced metric gij (x) , and thus the 
equality of Z with the partition function Z — J VXe^^^^ for the original energy ^ [or ([|)], if the parametrization is 
fixed likewise. In the sequel we shall measure the stiffness constants in units of temperature so that we can set T = 1 
everywhere. 

The advantage of the functional integral (^) is that the surface positions u°'{x) can be integrated out, leaving a 
purely intrinsic effective energy 

E^"" = ^ {Trlog [{D^ - D.X^W,] - £^ J d^x^gX^^h,,^ . (9) 

In the limit d — > oo, the fluctuations of the fields A'^ (x), hij (x) are frozen, and the energy is given by the saddle point 
approximation. For symmetry reasons, we may assume the saddle point to have constant diagonal fields hij — pdij, 
A*-' = A/i*-' — {X/p)6ij, so that the energy reads 

E - ^A..A.V{ [/ .og(t' + A.') - 2^a] + , (10, 

to be extremized in p and A. Here Aa;^Aa;^ is the base area of the surface. Performing the integral with a momentum 
space cutoff A of the order of the inverse molecular size ruv we find 



A <A W? '"^ + ^) = i t(A^ + A) log (A^ + A) - A^ - A log A] 



|fc|<A 

= i- [A2 (log A2 - 1) + (log A2 + 1) a - Alog a] + 0(1/A2), (11) 

and the brackets in ( |lO| ) become 

^ [2A2 (logA^ - 1) + (logA^ + 1) A - Alog A] - 2^ A. (12) 

The first constant term can be absorbed into the measure of the functional integral. The logarithmic divergence 
multiplying A may be removed by introducing a renormalized stiffness 



d - 2 / A2 



Ac^Ato-^^ (log— + 1) (13) 



where \x is some mass scale, on which k depends [k — ^(/i)]. For p equal to the inverse of the molecular size ruv, we 
have K(ruv) = '^o- With the help of k, the bracket in ( |iq ) becomes 
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1 A2 K 

Introducing further a /i- independent mass scale (the so-called dimensionally transmuted coupling constant) 

X ^ ^^2g-[2/(d-2)]47rK(M) + l^ 

expression ( p^ ) takes a fi- and cutoff-independent form 



and the energy ( p^ ) can be rewritten as 

d-2 



E = 



' f ^ ^ ^ 
Jo - -;- + K- 
47r p 



(14) 



(15) 



(16) 



(17) 



A multiphcative renormahzation factor = kq/k has been absorbed into A and A, so that all quantities are now 
finite for A oo. Extremizing ( p^ ) in p yields /o — A/47r and thus A = A, where /o — A/47r = 0. Extremizing ( p7| ) in 
A yields k./ p = l/47r, so that the extremal energy (p^) is 



d-2 ^ 1 A 2 ^ 
2 47r 



(18) 



From the trace of the logarithm in (|9|) we see that A sets a mass scale for the correlation function {DiU°- DjU ) which 

falls off like e^'^'^ for large \x\, showing that 1/VA plays the role of the persistence length (Q). Thus the d ^ oo 
model possesses precisely the properties which we derived for a real membrane in three dimensions by a one-loop 
approximation. 

The question of corrections to the one-loop approximation is therefore equivalent to the problem of lowering d down 
to the physical dimension 3. Since we are unable to treat the model for finite d, let us gain insight into its possible 
properties by comparing it to an analogous very similar model, which is also exactly solvable for a parameter N oo, 
which plays the same role as d — 2 in the above discussion. This is the O(A^) nonlinear cr- model for which the case 
N = 2 is well-known to have quite different properties from those derived from the N oo -limit. This knowledge 
will shed some light upon the behavior of surfaces for small d — 2. 



III. RELEVANT PROPERTIES OF 0(Af)-SYMMETRIC NONLINEAR a-MODELS 

The 0(A^)-symmetric nonlinear cr-model consists of a fluctuating field of unit vectors with N components 

na{x) (a = 1,...,7V), nl{x)^l, (19) 
in a two-dimensional x-space, with a partition function given by the functional integral 

d^-^n 



s 



N 



exp 



(20) 



where Si^ = 271^^^ /T{N/2) is the surface of a sphere in N dimensions covered by the directional integral d^~^n°-. 
This ensures a unit integral / d'^~^n°- / Sn — 1- Just as the surface model for d — > oo, this model is exactly solvable 
in the limit N —> oo, where it is referred to as the spherical model, first solved by Berlin and Kac in 1952 Pq |. As in 
the energy (H), we introduce an auxiliary field A(a;) and rewrite (|2^) as 



PAexp <j -y / d^x 



{duaf + NX{nl-l) 



where the A-integrations run from —ioo to ioo. Now the na-integrals are Gaussian and can be done, leading to 
Z„ = J mexp|-y Ti\og{-d'^ + X)+KQjd^xX{x) |. 



(21) 



(22) 
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In the hmit iV — > oo, the A-fluctuations are frozen at the saddle point, which lies at a constant X{x) 
energy 



N 1 , 
E„ = —Ax^Ax^ 

Regularizing the integral as before, this becomes 



A, and has an 



(23) 



E = ^Ax'Ax' p/o, /o ^ (log J - 1 



(24) 



where X = jj? exp{—A7T K{fi) / N} is the dimensionally transmuted coupling constant of the tr-model. It gives a nonzero 

length scale 1/Vx to the fluctuations of the vector field which have a correlation function c>c e~'^'^. 

The similarity of the two models is quite obvious. The advantage of the latter is that it is much better understood. 
Since the work of Kosterlitz and Thouless in 1973 [|l6| it is known that N ^ oo properties are found in the model 
only down to = 3 (classical Heisenberg model). For N = 2, the model possesses at large kq > Kq = 2/tt an extra 
phase in which the correlation functions of the vector field Ua have a long range, falling off algebraically like l/|x|™"*'' 
rather than exponentially. This is most easily seen by parametrizing the two-component vector field Na in terms of 
an azimuthal angle 9 as {Na) = (cos 6*, sin 0), and rewriting (pG) as a functional integral 



de{x) 



2tt 



exp <j / d'x[d,9{x) - 2ttS,{x; L)]'' 



(25) 



where Si{x; L) = J d£^j5^^\x — x{£^)) is the (5-function on a line L described by a;(^), pointing orthogonal to the line 
elements. The sum runs over a grand-canonical ensemble of lines L. The sum is necessary to preserve the cyclic 
invariance of the energy in ( pO| ) under cyclic replacements 0(.t) 9(x) + 2'k5(x\ S), where 5(a;; 5) = / d^^J^^^ (x— a;(^))) 
is the (5-function on the surface S described by x{^) jT^. The energy in ( p5| ) is a direct consequence of the identity 
Q^^^e(x) ^ e,''^(^)[d,9{x) - 2n5^{x]L)]. 

We may now introduce an auxiliary field bi{x) and rewrite Z„ as 



En 

L X 



d9{x) 

271 



dbi 



exp < Kq 



d^x 



bUx) 



ib,{x)[d^9{x)-2TTS,{x;L)] 



(26) 



The integrals over 9{x) can now be performed, which gives the condition that the field bi must be divergenceless [ p7| . 
This may be enforced by setting bi — eijdju, and the partition function ( |2^ ) becomes 



E 

{Xi,Xf} X 



n 



du{x) 



exp < Ko 



d^x 



[diu{x)Y — 2'Kiu{x)n{x) 



where 



n{x) 



,d^Sj{x;L) = S^^\x - x,) - S^'^\x - Xf) 



(27) 



(28) 



is the density of initial and final end points of the line L, which lie at Xi and x/, respectively, with positive and negative 
signs. These are the positions of vortices and antivortices in the original field configurations 9{x). By integrating out 
the u-ficld in (|27|), we obtain 



where 



Za = ^ exp 

{Xi,Xf} 



G{x - x') = 



— ^^TT^ / d^x I d^x n{x)G{x — x)n{x) 



d'^k e^^^ 



(27r' 



1 1 \x-x'\ riR 
2 i.2 =-^log^— +log — 



''0 



(29) 



(30) 



'0 



is the correlation function of the field u{x) containing some finite length scale rp and an infrared cutoff rm which 
goes to infinity. However, this infinity drops out in the partition function ( p9| ) since J d'^x n{x) = due to the equal 
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number of initial and final endpoints Xi,Xf in the sum over all initial and final points Y,xi.xf of the lines L in the 
grand-canonical ensemble. 

Let us derive the partition function ( p9| ) in a slightly different, more phenomenological way which will be useful to 
understand the smooth-membrane problem. Suppose we are given the partition function ( p7|) without vortices, where 
the energy is simply 

^lin Sxid^uf. (31) 

The equation of motion 

d^u{x) = (32) 
has two types of rotationally invariant solutions. The trivial u{x) = 0, and the solution 

1 r 

Uwh(a;) = T^Tj-log— , r > ro, r=\x\, (33) 

which satisfies — 9^Uwh = 0. The amplitude A is arbitrary. If we interpret u{x) as a vertical displacement of a flat 
surface in three dimensions, this looks like a "wormhole" in a curved planar space, as pictured in Fig. ^. 

u{x) 




This solution exists for any rg. Let us now imagine the surface to consist of molecules of size ryv- Then the smallest 
possible ro is ruv- For this value, however, the throat of the wormhole can be closed by a small cap from above or 
below, thus becoming a spike- like surface shown in Fig. ^to be denoted by Usp{x)^ or an antispike-like surface pointing 
downwards. 




FIG. 4. Spike-like solution of differential equation d'^u{x) = 0, obtained from the lower branch of the "wormhole" solution 
in Fig. ^ after capping it at the molecular scale ruv. 

The linearized energy ( ^2|) of these surfaces is large so that only a few spikes and antispikes are present in thermal 
equilibrium. If a single spike-antispike pair with opposite orientation centered at different places x and y is inserted 
into the field energy (|l|) , we obtain precisely the interaction energy of a vortex pair in the partition function ( ^ ) if 
we set A = \. The size of A is a consequence of the vortex quantization in the original energy in (|25|). A spike and 
an antispike can combine to a dipole, and two of these to a quadrupole, as shown in Fig. |^. In principle, there can 
also be a crystal-like array of spikes and antispikes, also shown in Fig. |^. This requires the presence off additional 
higher-gradient terms in the energy ( ^ which would modify the short-distance properties appropriately ||l^ . 
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The partition function (g9|) describes a neutral Coulomb gas. This is known to have a pair unbinding transition 
which is of infinite order at low temperatures, the Kosterlitz-Thouless transition. The transition point is easily 
calculated if only very few charges are present. At low temperatures, a single pair has an average square distance 

r°° 1 
^ ' Jr„ A-ko2tt ^ ' 

As the stiffness constant kq falls below Kq = 2/ti, the square distance diverges, and the electron gas enters a plasma 
phase in which the field u{x) acquires a finite range, the Debye screening length rocbyc; which is precisely the inverse 
of V\ in Eq. (H). 




FIG. 5. Pair of spikes, quadrupole-, and crystal-like multipole solutions, obtained from a superposition of spike and antispike 
solutions of the differential equation d^u{x) = in Fig. ^. 

This result implies that the average value of A(a:) which is calculated with the functional integral (p2|) as 



/ice f 
VX\{x) exp<^ - 
-ioQ L 



N 

y 



Tr log (-92 + A) + Ko / (fx \{x) 



(35) 



and which has the value (A(a;)) = A 7^ for small stiffness, vanishes for = 2 and kq > '*o- The fluctuations of \{x) 
become so strong that the nonzero saddle point at A = A looses control over the functional integral. The A(a;)-field 
will acquire more and more zeros until, at a small enough N and large enough stiffness, these zeros proliferate making 
(A(x)) = 0. Up to now, nobody has been able to derive this result from the functional integral (p5[), but the above 
Coulomb gas argument proves that it must be true. 

For a dilute gas of vortices, i.e., for a high activation energy and a small fugacity z = e~^''l^ of a vortex, the 
sum over all vortices and antivortices in the partition function ( p7| ) can be restricted to a single vortex or antivortex 
at each point, in which case the partition function goes over into 



n 



du{x) 



exp 



d^X 



1 z 

— [diu{x)]'^ cos[27rKou(a;) 

2 Ko 



(36) 



This is the partition function of the sine-Gordon model, which is thus equivalent to the 0(2) nonlinear cr-model. For 
small stiffness kq, we may expand cos[27rKoM(2:)] ~ 1 — 'in'^ KqU^ (x) / 2 to see that the fluctuations have a a finite range. 
For large stiffness kq > 2/tt, however, the cosine oscillates so rapidly that it no longer contributes to the functional 
integral. In this regime, the system has only long-range fluctuations [18|. 
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IV. CONSEQUENCES FOR MEMBRANES 



What are the lessons of all this for membrane fluctuations? Suppose we restrict the field fluctuations in (^) to 
diagonal fields hij{x) — p{x)Sij and A'-' = X{x)h''^ (x) with the constant extremal p = At:k. Then we remain with 
a functional integral over X{x) with a Boltzmann factor where E'°^ is the obvious generalization of the 

effective energy (|9|) to a functional of a space-dependent A(a;) which has precisely the same form as functional the 
exponent in But such a functional integral has just been shown to yield for a large enough stiffness k a vanishing 
average (A(x)), implying that the fluctuations of the derivatives 9^?/° have an infinite range characteristic for a smooth 
surface. It is certainly conceivable that the neglected fluctuations, nondiagonal in A*-' and arbitrary local in hij (x) do 
not change this result. 

How can we see whether this is true? Consider a linearized version of the energy m) 



(37) 



and let us go through the same argument as in the discussion of the linearized energy (^). The equation of motion 

d\{x) = (38) 

has three types of solutions. First, there is the trivial one ^(a:;) = 0. Second, there are spike-like solutions Usp{x) 
arising from capping a wormhole solution ( |3^ ) of the previous system, as shown in Fig. |[ Third, and most importantly 
for us, there are spike-like solutions of the form 



Usp2[x) = ±A— — 

Z OTT 



i log log ■ 

roe ? 



roe 



1/6 



(39) 



which satisfy — 9^Usp2(a;) — Usp{x) and 9^Usp2(a^) = ±v4(5^^-' (x). With a cap of molecular size they have the form 
shown in Fig. |^. There are also dipole, quadrupole, and multipole solutions displayed in Fig. 0. 




FIG. 6. Spike-like solution of equation d*'u{x) = 0, with singularity capped at the molecular scale ruv. 

The distances of the spikes in these arrays are of the order of the length scale r,n ~ 1 /cm supplied by the extrema 
in Fig. |. 



FIG. 7. Pair of spikes, quadrupole, and multipole solutions of differential equation d'^u{x) — 0, with singularities capped at 
the molecular scale ruv- 



The important point is now that there must be a preferred amplitude A of the spike- like solution (|39|), fixed by the 
interplay between and the molecular size ryv (as in the Monte-Carlo calculation in Ref. Q). If this is the case, 
we can now insert a superposition of these solutions into the energy (pT\) and find a partition function 



J2 

{Xi,Xf} 



-A' 



Sx I d^x' n[x)Gi{x — x')n{x') 



(40) 



where n{x) is a sum of (5-functions at the positions of the spikes (with a minus sign for antispikes), and 

(Pk e''^!^-^') 



Ga(x — x) = 



(27r) 



fc4 



— Usp2{x — X) + Ci{x — x) + C2 



(41) 



is the correlation function of the field u{x) containing two infinite constants from the infrared divergence of the 
integral. The infinity C2 drops out for neutral spike-antispike systems with J (Pxn{x) = 0, the infinity ci drops out 
under the condition of dipole neutrality J d'^xxn{x) — 0. The latter follows directly from 



d X d x' n{x) {x — x) n{x) — 2 I d xx n{x) I d x n{x') 



d xxn{x) I d X X n{x') 



= -2 / d^xxn{x) I d^x' x' n{x') = Q. 



(42) 



Note that the simpler spike- like solutions Usp{x) in Fig. ^ do not have a long-range attraction between them. 
The partition function (EO) is obviously equivalent to the local one 



or, in analogy with ( 

Zn 



E n 



du{x) 



exp < Kq d X 



[d'^u{x)]'^ - Aiu{x)n{x) 



(43) 



to 



n 



du{x) 



exp 



-Kq 



d^X 



1 Z 

-[9^M(a;)]^ cos\Akou{x)] 

2 Ko 



(44) 



This system is the dual transform of a two-dimensional crystal with defects whose phase properties have been thor- 
oughly studied |Q. As it stands, the partition function ( p^ ) has a single first-order transition j2l] from a spiky to 
a smooth surface. If the short-range properties of the model are slightly modified at a length scale Z, the first-order 
transition splits into a sequence of two Kosterlitz-Thouless-like transitions, in which the spikes first form a gas of 
dipole pairs, and then a gas of quadrupoles (see the figure on p. 1303 of Ref. ||20[|). 

Upon further modifying the short-range properties, the quadrupoles may combine to a liquid of spikes and antispikes, 



which may be the structures seen in the experiments in Fig. ^J. The liquid could also freeze to a solid. Future 
experiments should look for quadruplets and dipoles to confirm the correctness of these ideas. 

It is essential that nonlinearities in the curvature energy fix an optimal size of the amplitude A in (40). For a 
purely quadratic surface energy (|l|) such an optimal A cannot exist since the energy of a spike-like solution ( |39| ) can 
be reduced continuously to zero by letting A go to zero. 
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V. POSSIBLE CONSEQUENCES FOR STRING THEORY 



Such nonlinearities may also resolve an an old riddle in the string theory of permanent quark confinement. Quarks 
and antiquarks are held together at all distances by color-electric flux tubes formed from nonabelian gauge fields of 
quantum chromodynamics. In a euclidean spacetime with d = 4, these tubes form fluctuating surfaces. In string 
theory, the behavior of such surfaces is investigated in an idealized form, considering them as infinitely thin objects 
possessing initially only surface tension {Namhu-Goto strings). But soon it was found that these surfaces would be 
unstable against the formation of infinitely thin protrusions, which do not possess any surface energy but a large 
configurational entropy [plumber's nightmare). This instability would ruin quark confinement. To avoid this, a 
curvature stiffness energy was added to the surface tension p3t . The resulting stiff-string action was equal to the 
energy (^. However, this did not seem to lead to the desired stabilization since, in an infinite number of dimensions, 
there exists a finite persistence length beyond which quark would no longer be confined. It remained unclear whether 
this phenomenon would persist down to d = 4. The present discussion gives rise to the hope that the persistence length 
can become infinite after all. Since d = 4 in (||) corresponds to iV = 2 in the nonlinear cr-model (p2[), which possesses 
a smooth phase with long-range correlations, the A-fluctuations in the action (||) may be violent enough to produce 
so many zeros in X{x) that (A(a;)) — and the persistence length in infinite, implying permanent quark confinement. 
But we also have learned that this may occur only with support from higher nonlinear terms in the curvature energy. 
A color-electric flux tube certainly has such additional terms, although it differs from the membranes discussed in this 
l^a^ er in an important point: Its curvature stiffness is negative, as was recently shown by the author and Chevyakov 

Paper s by the author in the list of references can be read on the internet, for instance (pubsl#181) under the www 
address trttp : / / www. physik . fu-berlin . de/ ~ kleinertkleiner _rel ^ 181 . 
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